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INTRODUCTION 

Theories based on Supersymmetry (SUSY) are widely 
considered as the theoretically most appealing extension 
of the Standard Model (SM). The Minimal Supersym- 
metric Standard Model (MSSM) predicts the existence 
of scalar partners /lj/a to each SM chiral fermion, and 
spin-1/2 partners to the gauge bosons and to the scalar 
Higgs bosons, where two Higgs doublets are present in 
the MSSM. So far, the direct search of SUSY particles 
at present colliders has not been successful. One can 
only set lower bounds of o(lOOGeV) on their masses, 
see Ref. (1). An alternative way to probe SUSY is to 
search for the virtual effects of the additional particles. 
The experimental precision has to be matched with high 
precision theoretical predictions for the various precision 
observables. The most prominent role in this respect is 
played by the p-parameter, see Ref. (2): 
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The radiative corrections to the vector boson self-energies 
at zero momentum transfer, Ez.iv(O), constitute the 
leading, process independent corrections to many elec- 
troweak precision observables, such as the W boson mass, 
Mw, where 



mw^Mw/2 4/(4- 4) Ap, 
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with 4 



1 - 4 = M^/Ml. Within the MSSM the 
corrections have so far been restricted to o(aav), see 
Ref. (3). In order to match the accuracy obtained in the 
SM and the (prospective) experimental uncertainties, the 
leading two-loop corrections to Ap at O {G^mf) are desir- 
able. These corrections are the lowest order contributions 
involving non-SM particles with a mf dependence. In the 
limit of a large SUSY scale, Msusy > Mz, the SUSY 
particles decouple from the observables, leaving the two 



Higgs doublets of the MSSM active. First results for this 
contribution have recently been calculated in Refs. (4, 5). 

Feynman-diagrammatic two-loop calculations involve 
a large number of diagrams. In the MSSM the additional 
problem of a prohferation of scales is apparent. There- 
fore, in order to perform the calculation as presented in 
Ref. (4, 5), the use of computer algebra programs is in- 
evitable. In particular, we made use of the amphtude gen- 
erator FeynArts, see Ref. (6), where the MSSM model file 
has recently been completed. The reduction of the am- 
plitudes to scalar integrals have been performed with the 
program TwoCalc, see Ref. (7). 



TECHNIQUES FOR TWO-LOOP 
CALCULATIONS 

In order to calculate the {Gpmf) corrections to the 
p-parameter, the diagrams in Fig. 1 have to be evaluated. 





FIGURE 1. Generic diagrams for o{Gpmj) corrections to Ap. 
(V=Z,W, f = t,b, = h,H,A,H^,G,G^) 

The diagrams and the corresponding amplitudes have 
been generated with the help of the Mathematica program 
FeynArts and the recently accomplished MSSM model 
file. This model file contains all relevant information 



about the MSSM particles and vertices. (Only MSSM 
counter term vertices are missing at present.) FeynArts 
has been checked in several ways to insure its rehabil- 
ity. Apart from many self-energy calculations, whole 
processes like e+e^ tt, e^e~ H^H^ , qq tt and 
e+e" X^'^x'^'^ have been calculated by hand and with 
FeynArts. Perfect agreement has been found for all pro- 
cesses, see Refs. (6, 8). 

How easy the evaluation of the amplitudes has be- 
come with FeynArts is demonstrated in the following se- 
quence from a Mathematica session where the Z boson 
self-energy amplitude, corresponding to Fig. 1, has is ob- 
tained. (A detailed guide can be found in Ref. (6).) 

(start of Mathematica session) 

> <<FeYnArts .m; 

{-^ loading FeynArts into Mathematica) 

> se2 = CreateTopologies [ 2 , 1->1, 

ExcludeTopologies->Internal ] ; 
{-^ two-loop topologies are created) 

> V2V2 = InsertFields [se2, 

V[2] -> V[2], 
Model->"MSSM", 

ExcludeParticles -> { . . . } ] ; 
(^ fields are inserted into the topologies, incoming 
field V [ 2 ] = Z and outgoing field V [ 2 ] are specified, 
the model is chosen to be the MSSM) 

> Paint [V2V2] ; 

(^ the diagrams are painted (optional), see Fig. 1) 

> V2V2A = CreateFeynAmp [V2V2] ; 

(^ Feynman diagrams are converted into amphtudes) 

> V2V2A » V2V2.amp; 

(— > result for the two-loop amplitude of the Z self- 
energy is saved in the file V2V2 . amp) 
(end of Mathematica session) 

The further evaluation of the amplitudes has been 
performed with the Mathematica program TwoCalc, see 
Ref. (7). It performs the reduction of the self-energy am- 
plitude at the two-loop level to a basic set of scalar inte- 
grals. The appUcation of TwoCalc is demonstrated in the 
following sequence from a Mathematica session where 
the Z boson self-energy amphtude is processed. 

(start of Mathematica session) 

> <<TwoCalc.m; 

(^ loading TwoCalc into Mathematica) 

> amp =<<V2V2 . amp; 

(— > the amplitude obtained with FeynArts is loaded) 

> SetOptions [ TwoLoop, 

CollectFunction -> 01; 



(-^ options are set. CollectFunction allows to 
choose between DREG and DRED) 

> res = TwoLoopSum [ amp, 

SelfEnergyPart -> 1]; 
(— > the amplitude is reduced to scalar integrals) 

> res = Collect [res, 

{AO [_] , T [ ] } , Simplify] ; 

(^ the result is simplified (optional)) 

> res >> V2V2.res; 

(— > the result for the Z boson self-energy amplitude in 

terms of scalar integrals is saved in V2V2 . res) 
(end of Mathematica session) 

The analytically obtained result now consists of scalar 
integrals at the one- and at the two-loop level. In the 
case presented here, since the external momentum is set 
to zero, these are the functions Aq and T134 at the one- 
and two-loop level, respectively, see Refs. (9, 10). 

The pure two-loop diagrams have to be supplemented 
with the corresponding one-loop diagrams with counter 
term insertion. In order to obtain the 0{Gpmj) correc- 
tions an expansion of the amplitudes up to 0{mf/M^) 
had to be performed. All remaining Mw = Mzcw have 
been set to zero. Furthermore we had to apply the MSSM 
sum rules for Higgs boson masses, especially implying 
for the case Mz.Mw —> that the lightest MSSM Higgs 
boson has the mass m/j = at tree level. After adding up 
all contribution, the expansion of the result in terms of 
5 = (4 — Z))/2 leads to a finite result in the limit 5 — > 0. 
All relevant details can be found in Refs. (4, 5). 



RESULTS FOR Ap AND Mw 

After extracting the prefactor of mf /M^ and setting 
Mw to zero, besides sp = tan p/ ■\/ 1 + tan^ p only two 
mass scales remain: the top quark mass, and the mass 
of the CV-odd Higgs boson. Ma. The result for Ap can 
be conveniently expressed in terms of a = mj/M^: 
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with A = 3— 13a + llfl^. As a consistency check, in the 
limit of Ma ^°°,a^O,wc obtain 
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This is the SM result with Mh — > 0. It shows that the 
MSSM decouples to the SM limit (also at the two-loop 
level) when the new scales, here Ma, is made large. 
In the limit of large tan P (i.e. (1 — ip) <C 1) one obtains 
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Thus for large tanP the SM limit with Mh ^ is reached. 

The numerical result obtained for the O {Gpirif) cor- 
rections to Ap vary around —2 x 10^^ and are of the 
same size (but with opposite sign) as the leading o (aa^) 
MSSM corrections originating from the scalar top and 
bottom sector, see Ref. (3), if the common soft SUSY- 
breaking scale is chosen to be Msusy ~ 500 GeV. For 
a larger value, Msusy ~ 1000 GeV, the new 0{Gpmf} 
corrections are about three times larger than the o(aaj ) 
contributions. Furthermore it is well known that the 
{G}mf) SM result with MfJ^ = underestimates the 
result with M™ 7^ by one order of magnitude. One 
can expect (see Ref. (4)) a similar effect in the MSSM 
once higher order corrections to the Higgs boson sector 
are properly taken into account, which can enhance nth 
up to mft < 130 GeV, see Ref. (11). 

With the help of eq. (2) the shift in the W boson mass 
can be evaluated. In Fig. 2 the induced shift for Mw is 
shown as a function of Ma for tan P = 3 , 40. The effect on 
5Mw varies between — l.SMeV and — 2MeV for small 
tanP and is almost constant, 8Mw = — 1.25 Me V, for 
tan P = 40. The constant behavior can be explained by the 
analytical decoupling of tanP in eq. (3) when tanP 1, 
see eq. (5). 

CONCLUSIONS 

We presented the calculation of the SUSY contribu- 
tions of o{Gj^mf) to the p-parameter. In order to ob- 
tain these two-loop corrections, the computer algebra pro- 
grams FeynArts and TwoCalc have been used. Examples 
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FIGURE 2. The induced shift in My/ is shown as a function of 
Ma for tanp = 3,40. 



of the handling of these programs were given. A com- 
pact analytical result in terms of mj/M^ has been de- 
rived. The induced shift in Mw varies between —1 MeV 
and —2 MeV, depending on Ma and tanp. 
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